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Abstract 

In this letter we study a novel effect of a hidden sector coupling to the standard 
model Higgs boson: an enhancement of the Higgs pair production cross section near 
threshold due to bound state effects. After summing the ladder contributions of 
the hidden sector to the effective ggHH coupling, we find the amplitude for gluon- 
gluon scattering via a Higgs loop. We relate this amplitude to the double Higgs 
production cross section via the optical theorem. We find that enhancements of 
the 0(100) for the partonic cross section near the threshold region can be obtained 
for a hidden sector strongly coupled to the Higgs boson. The corresponding cross 
section at the LHC can be as large as O(10) times the SM result for extreme values 
of the coupling. The detection of such an effect could in principle lead to important 
information about the hidden sector. 
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1 Introduction 



The idea of a hidden sector, that is, a sector that is singlet under the symmetries of the 
standard model and interacts with the rest of the particles only through the Higgs boson 
(in addition to gravity), sometimes called "Higgs portal", was probably first introduced 
by Veltman and Yndurain [TJ with the purpose of parametrizing the limit of large Higgs 
mass in the standard model. 

A more recent motivation of a hidden sector comes from new viable dark matter 
models. The existence of dark matter inferred from several different observations signals 
the incompleteness of the standard model of the electroweak interactions [2]. Perhaps the 
best motivated candidate for dark matter are neutralinos arising from supersymmetric 
extensions of the standard model but there are other contenders from different extensions 
with varying degrees of theoretical motivations [3j. 

The simplest extension of the standard model with a natural dark matter candidate 
is the addition of a scalar singlet [I]. A more recent role of singlet scalars in dark matter 
phenomenology appears when one tries to explain the excess of positrons and electrons 
seen in some experiments as an effect of dark matter annihilation [3]. In this case, the 
annihilation cross section must be a factor of O(10 3 ) larger than the usual cross section 
determined by the observed dark matter abundance in the case of thermal relics. The 
existence of new light singlet scalars coupling to the dark matter candidate can generate 
an enhancement in the annihilation cross section of dark matter particles, known as 
the Sommerfeld enhancement [5J. Several studies have been performed considering this 
possibility [6]. In particular, the coupling of this scalar to the Higgs boson can lead to 
sizeable effects in the direct detection of dark matter, as recently emphasized in [TJ. 

In this letter we want to study a novel effect of the coupling of the Higgs boson with 
another scalar field: the enhancement of the double Higgs production at the LHC near 
threshold due to the formation of Higgs-Higgs bound states. 

Double Higgs production is an important process to test the structure of the Higgs 
boson potential and possible new physics beyond the Standard Model [8]. A possible 
enhancement effect may be of importance since the gluon distribution functions are largest 
near the threshold region. This is akin to the well known enhancement of top quark pair 
production near threshold due to the gluon contribution [9]. Contribution of bound states 
to the production of supersymmetric particles at the LHC has also been recently discussed 
in [10]. The enhancement of the double Higgs production may have observable effects 
and therefore can in principle probe the hidden sector. 

2 Enhancement 

The formation of a nonrelativistic Higgs-Higgs bound state, sometimes referred to as Hig- 
gsonium [TT] or Higgsium [12] would result in an increase of the double Higgs production 
cross section near threshold. The possibility of formation of Higgsium can arise within 
the standard model due to the Higgs boson self-interactions, but only in the case of a very 
heavy Higgs, actually above the unitarity bound [13J. Since we will be interested in an 
intermediate Higgs boson mass, we will only consider the contribution of the interaction 



of the Higgs boson with the hidden sector to this process. 

Double Higgs production within the standard model via the dominant gluon fusion 
process at hadron colliders was computed by Plehn et al. [14J. They use a form-factor for 
an effective gghh coupling that we will denote by T^. This effective coupling depends 
dominantly on the total momentum entering the vertex. 

We add to the standard model lagrangian a new singlet real scalar with a coupling 
to the Higgs boson h: 

£ = £sm + - (<9 M 0) 2 - ^mj0 2 + guhhcj) (1) 

where m^ is the mass of the new scalar, g is a dimensionless coupling constant and u 
is an energy scale from this new sector. This lagrangian should be thought of as a toy 
model used to study a general phenomenom. For the moment we will assume that it is 
already in its diagonal basis for the scalar sector. 

Unfortunately there is no full treatment of bound states in quantum field theory and 
some approximations must be used, such as the ladder approximation to the Bethe- 
Salpeter equation [15]. We will follow the standard framework developed for computing 
threshold effects due to bound states [9] and use the optical theorem to write: 

a(gg ->• hh) = -ImM(gg -> gg), (2) 

S hh 

where Ai(gg — > gg) is the amplitude for gluon-gluon scattering through a Higgs boson 

hh 

loop. Our goal will be to compute this amplitude in the nonrelativistic limit taking into 
account the contributions from the hidden sector. 

We first compute the corrections to the effective gghh vertex due to the hidden scalar 
exchange by summing all ladder contributions as shown in Figure (lit), which results in 
the integral equation 

= 1 + 14^ , 'ff 1 (3) 



(2tt) 4 (p-kf 



m 



S(q,k) 



(q/2 + k) 2 -m 2 h - im h T h (q/2 - k) 2 - m\ - im h T h 

where S(q,p) = T(q,p) /T^ ' (q,p) is the normalized coupling, q is the 2-particle bound 
state momentum and 2p is the relative momentum 0. A bound state would appear as a 
pole at q 2 = m 2 hh , with rrihh < 2mh- Near threshold we write q = (2rrih + E,0), where 
\E\ « rrth is the binding energy for a nonrelativistic bound state. 

In the nonrelativistic limit we keep only the instantaneous part of the <fi propagator. 
Likewise, since we are in the threshold region, we neglect the time dependence of the 
total vertex gghh and for consistency we do the same for the vertex inside the loop in 
the Bethe Salpeter equation (figure [T^,) . With these approximations one can perform the 
k integral to obtain: 

g 2 u 2 f d 3 k 1 1 



om h J (27TJ (p — k) 2 — mi k 2 /m h — E — %T h 



1 One should note that a including a hhtficj) coupling in eq.([T]) would result in an additional contribution 
with a loop of (j> held. This contribution would be suppressed by the usual g 2 /16n 2 factor. 
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Setting 



G(q,p) 



S(q,p) 



(5) 



fp/nrih — E — iT h 

and performing a Fourier transform in Eq.(j4]) we arrive at a Schroedinger-like equation, 
depending of the relative position r and center of mass coordinate r' placed at origin: 



E - iT h + V(r) 



G v (r,r =0,E) =5 3 {r) 



(6) 



with an Yukawa potential given by 

V(r) 



g 2 u 2 e -m^r 

8m\ r 



(7) 




b) 



hh 




Figure 1: (a) Bethe Salpeter equation for Ai(gg — > hh) with hidden sector corrections 
in the ladder apro 

and M.{gg —> hh). 



in the ladder aproach (b) Relation between the total scattering amplitude M.(gg — > 

hh 



Now we can find the amplitude M.{gg —> gg) by solving the integral equation depicted 

hh 

in Figure (Jib): 



M(gg->gg) = T<®(E) 

hh 

d 4 k 



T(q,k) 



(2tt) 4 (q/2 + kf -m 2 h - im h T h (q/2 - k) 2 -m\- im h T h 
Following similar steps as in the previous calculation we obtain 

M(gg^ h gg) = {T^(E)) 2 G v (0, E). (9) 

Therefore the cross section with the contribution from the hidden sector is 

a(gg -> hh) = a (gg -)■ hh)R(E) (10) 
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with the enhancement factor R(E) is given by 



lmG (0,E) 

where E = y/s — Im^ is the center-of-mass energy of the Higgsium from threshold and 
G is the solution of the same Schroedinger equation (j6]) in the absence of a potential. 



3 Finding the enhancement factor 

The solution to eq.(|6]) can be obtained following the standard procedure of defining two 
independent solutions V\ and v 2 of the corresponding radial homogeneous equation in the 
relevant case of zero angular momentum, one regular at the origin and the other regular 
at infinity: 



Air rr 

where Vi^{r) is a solution of the equation: 



for < r < r' < 00 
for < r' < r < 00 



G(r,r',E + zT h ) = { ^^Cm , „ , . (12) 
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l_ + m , (E + iT h - V{r)) 



v(r) = (13) 



We use the method described in Strassler and Peskin [9] to numerically solve the 
equation above with the appropriate boundary conditions. We will use the functions v a 
and Vb with boundary conditions v a (r — > 0) = r and (r — > 0) = 1 to write 

vi(r) = v a (r) (14) 
v 2 (r) = v b (r) + Bv a (r). (15) 

Hence B = lim (—^t\) and we can formally calculate the imaginary part of the Green 
: 

ImGV(0, 0, E + iT h ) = -^ImB (16) 

47T 

In our numerical solution the boundary conditions for the Green function is of course 
calculated at a finite value of r, and it is very important to keep the Higgs boson width 
that guarantees its exponential decay. 

In the absence of a potential it is easy to find that 

Go (ry,E + ir,) = -^^l e — (17) 

47T Ar r 

where A = yJm(E + iT h ) and therefore 

1mG (0,0, E + iT h ) = -^ReA. (18) 

4"7T 



2 Note that the real part of this limit diverges. 
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4 Results and Conclusion 



We use the results of Plehn et al. [H] to evaluate the cross section without hidden sector 
corrections 0, which can be written as: 



a (gg -> hh) = -L \ M (gg -+ hh)\ 2 ^ (l - ^ . (19) 

Part of the kinematical factor is cancelled by the denominator of the enhancement 
factor of eq. tfTO]) resulting in a total cross section given by 




a(gg -> hh) = — \M(gg hh)\ z J ( I ) -== (20) 

The enhancement depends basically on two parameters: the strenght of the interaction 

2 2 

k = and its range determined by the mass of the exchanged boson. For = 

h 

one is back to the well known Coulomb case. As an illustration of this effect we will set 
the Higgs mass in 180 GeV (we checked that our results are weakly dependent on the 
Higgs mass in the intermediate range). 




V7- 2 m h (GeV) V7 - 2 m h (GeV) 

Figure 2: Double higgs production cross section as a function of energy above threshold 
for m<£ = 50 GeV. Left panel: Double higgs production cross section cross for k = 
0.01,0.1, and 0.3 (dashed, dot-dashed and dotted lines, respectively) compared with the 
uncorrected result (solid line). Right panel: Double higgs production cross section near 
threshold for k = 0.1,0.7, and 1 (dashed, dot-dashed and dotted lines, respectively) 
compared with the uncorrected result (solid line). Notice the logarithm scale in this case. 



In Figure (|2]) we show the total cross section a(gg — » hh) for different values of the 
coupling constant for = 50 GeV. For couplings k < 0.01 hardly any modification is 
observed. The difference below threshold is due to the finite width effects contained in 
the function G, which allows a nonzero cross section below threshold. We have checked 
that our numerical code reproduces the analytical result eq. (fT8|) to high accuracy. 

3 We will use the y/s » m t limit for illustration. 
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Table 1: Enhancement factor for double Higgs production cross section at the LHC for 
different values of the coupling k to a hidden sector for \fs = 7 and 14 TeV 

The enhancement increases with the coupling and for values k > 0.6 a dramatic dif- 
ference results from the presence of peaks due to Higgs-Higgs bound states. We find 
enhancement factors as large as a factor of (9(100) in this case, which is comparable to 
what has been found in the study of Sommerfeld enhancement for dark matter annihila- 
tion cross section in several models [5]. 

A couple of comments are in order. First, it is a textbook exercise [16] to implement a 
variational method to estimate the energy of the 1=0 bound state in a Yukawa potential 
In our case the energy is given by: 

E = - K m^ 3 (p-l)/4(p+l) 3 (21) 

where p is a solution of the equation Krrih/m^ = (p + l) 3 /p(p + 3). For k = 1, = 
180 GeV and = 50 GeV, the values used in our paper, we obtain E = —11 GeV 
(E = —1.8 GeV for k = 0.7). These are in good agreement with our numerical results 
shown in Figure(|2]) and corroborate that the bound states are nonrelativistic, since the 
binding energy is much smaller than the Higgs mass even for these large values of kappa. 
Secondly, the nonrelativistic approximation is only valid for a region of 30 GeV or so 
around threshold. The left plot of Figure (T5J) is only meant to illustrate that the corrections 
are small even at high energies. In our numerical results for the LHC cross sections below 
we switched-off the bound state corrections at an energy of 15 GeV above threshold. 

For lighter masses (O(10) GeV) the peaks appear for smaller couplings at lower en- 
ergies and occasionally more than one peak can be seen. For heavier masses, as expected, 
bound states are not formed but some enhancement effect is still obtained. 

In order to estimate the corresponding enhancement at the LHC, we have convoluted 
the parton level cross section with the gluon distribution function in the proton. For our 
estimates we used the Mathematica package of the leading order CTEQ5 [17] with fac- 
torization and renormalization scales set to q = 2rrih. In Table (II]) we present our results 
for \fs = 7 and 14 TeV. The convolution smooths out the large partonic enhancements 
but one can still find cross sections that are as large as 20 times the SM result for extreme 
values of the coupling. We find that the enhancement factors are not very sensitive to the 
CM energy, although the absolute values of the production cross section certainly are. 

At this point we should note that the simplest model with an additional singlet would 
not result in a large enhancement due to the mixing between the scalars, which requires 

n < where A is the Higgs boson self-coupling [18]. There are a few examples where 



J We thank G. Krein for pointing this out to us. 
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one could have a strong interaction of the hidden sector with the Higgs boson, such as 
considering a more evolved scalar sector or even interactions with hidden abelian gauge 
fields or a condensate of fermions. We take no position on what the primary motivation 
might exist for the hidden sector. Our goal is to consider the general effect of this hidden 
sector on double Higgs production near threshold. 

One may also be concerned with unitarity violation when large couplings are present. 
For example, at low energies one gets a s-channel contribution to the s-wave scattering 
amplitude that goes as ao = g 2 u 2 /327rm|; this is less than one for g 2 u 2 = 0{m\) and the 
examples discussed here. 

The prospects for detection of such an enhancement require a detailed study that is 
beyond the scope of this short letter. The signal would probably be an excess of 4 gauge 
boson events near the double Higgs threshold, where the SM background is small. It is 
amusing to entertain the idea that an investigation of the double Higgs cross section near 
threshold could lead to important information about the hidden sector. 
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